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摘 要 : 基于 行星 会 合 指数 运动 学 方程 ,推演 出 太阳 质心 到 太阳 系 质心 距离 变化 的 计算 公式 ,并 利 
用 该 公式 重建 太阳 质心 相对 于 太阳 系 质 心 的 距离 变化 时 间 序 列 , 作 为 太阳 轨道 运动 的 蔡 代 性 指 
标 , 使 用 聚合 经 验 模 态 分 解 (Ensemble empirical mode decomposition , EEMD ) 方 法 分 析 0~25.0 ka BP 
太阳 质心 到 太阳 系 质 心 的 距离 与 A"C 时 间 序 列 。 结 果 显 示 : 两 者 均 存 在 蕴 合 ~2300 a 周期 (Hall- 
statt cycle) 的 IMF 分 量 ; 互 相关 分 析 证 实 两 者 列 含 的 Hallstatt 周期 分 量 间 具有 相关 性 ,0~13.9 ka BP 
与 13.9~25.0 ka BP 年 代 相 关系 数 分 别 可 达 0.52、0.44, 并 且 太 阳 质 心 到 太阳 系 质心 距离 的 变化 超前 
于 太阳 活动 指标 A"C 的 变化 。 表 明太 阳 质 心 远离 太阳 系 质 心 时 ,大 气 "C 含 量 增加 ,对 应 于 太阳 活 
动 较 弱 的 时 期 ;太阳 质心 靠近 太阳 系 质心 时 ,大 气 "C 含量 减 小 ,对 应 于 太阳 活动 较 强 的 时 期 。 这 一 
认识 对 于 进一步 探究 太阳 活动 Hallstatt 周期 的 成 因 机 制 提 供 了 新 的 依据 和 思路 。 
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太阳 活动 的 强 弱 直接 或 间接 影响 地 球 系统 动 
力学 过 程 ,对 许多 地 球 物理 现象 影响 深远 ""。 因 此 ， 
对 于 太阳 活动 变化 规律 的 研究 一 直 受 到 学 者 们 的 
关注 和 重视 ”1。 

“C 生产 率 和 太阳 活动 强度 变化 有 关 “ ,其 作为 
太阳 活动 棕 代 性 指标 被 广泛 应 用 于 太阳 活动 变化 
规律 的 研究 中 。Suess“ 利 用"*C 检 测 出 太阳 活动 周 
期 2400 a; Damon 和 Sonett "检测 出 周期 2241 a; 
Vasiliev 和 Dergachev "检测 出 周期 2400 a; Fra IR 
等 中 对 Solanki 等 ”利用 "C 重建 的 太阳 黑子 数据 进 
行 Lomb-Scargle 谱 分 析 , 检 测 出 周期 2270 a。 但 是 ， 
关于 这 个 显著 的 ~2300 a 太阳 活动 周期 (Hallstatt cy- 
cle) 产 生 的 原因 ,还 缺乏 相对 一 致 的 意见 。 

对 于 太阳 活动 变化 机 制 目前 主要 有 两 种 解释 ， 


验 模 态 分 解 (EEMD ) 


系统 驱动 理论 。 太 阳 系 行星 系统 驱动 理论 ,又 可 分 
为 两 种 ,一 是 行星 系统 对 太阳 的 引 潮 力作 用 '"， 
二 是 行星 系统 作用 下 太阳 绕 太 阳 系 质心 的 轨道 运 
动 "1, 后 者 近年 来 逐渐 成 为 学 者 关注 的 焦点 。 
Charvatova ' 首次 提出 2100~2500 a 的 振荡 可 能 
是 由 于 太阳 绕 转 运动 ,其 利用 一 个 简化 模型 模拟 的 
太阳 绕 转 运动 ,具有 平均 2402.2 a 周期 ,相当 于 木 
星 - 日 心 - 地 心 排 成 一 列 (9.885x243), 在 每 一 个 周 
期 中 ,都 具有 370 a 的 标准 三 叶 形 循环 模式 ,于 
8000 a 的 "CC 代理 记录 中 对 应 太阳 活动 较 高 时 期 。 
Scafetta 等 1 计算 行星 系 质心 相对 于 太阳 运动 轨道 
的 瞬时 偏心 矢量 变化, 发现 具有 2318 a 周期 ,与 利用 
谐 波 曲线 拟 合 "C 去 趋势 分 析 后 的 残 差 记 录 对 比 , 图 
像 基 本 对 应 ,相位 差 5/2, 由 此 认为 在 Hallstatt 周 期 


一 种 是 太阳 发 电机 理论 ; 另 一 种 是 太阳 系 行星 
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的 形状 演化 为 更 圆 的 形状 时 ,产生 的 放射 性 核 粒 子 
会 更 多 ,从 而 调节 进入 太阳 系 及 地 球 的 宇宙 射线 通 
量 。 刘 复 刚 等 站 基于 修正 系数 (g=1.11518) 计 算 行 
星 会 合 指数 图 像 与 8000 a H “C 曲线 对 比 ,探讨 两 者 
之 间 的 关系 ,认为 行星 会 合 指数 极 大 值 包 络 线 中 的 
合 值 表明 此 时 太阳 所 具有 的 势能 较 小 , 癌 行 星际 空 
间 释 放 的 能 量 相对 减少 ,从 而 使 进入 地 球 大 气 的 宇 
宙 射 线 和 能 量 粒子 沉降 增加 ,”C 含量 增 大 。 

前 人 研究 已 取得 一 些 成 果 , 但 是 在 探究 太阳 轨 
道 运动 与 "C 关 系 时 ,采取 简单 图 像 对 比 ,或 预先 设 
定 函数 关系 式 进行 曲线 拟 合 ,于 分 析 结 果 注 入 了 人 
为 因素 ; 且 研 究 年 代 大 多 集中 在 全 新 世 ,而 放射 性 
WIER (INTCAL1I3 ™ ) 具 有 10 a 分 辩 率 的 A“"C 记录 
已 拓展 到 25.0 ka BP. 

2013 年 , 刘 复 刚 和 王建 六 提出 行星 会 合 指数 天 
研究 行星 系统 运动 规律 ,本文 基于 这 个 指数 ,推算 
出 太阳 质心 到 太阳 系 质心 的 距离 公式 ,并 计算 其 随 
时 间 变 化 的 序列 ;利用 EEMD 方法 对 0~25.0 ka BP 
太阳 质心 到 太阳 系 质心 的 距离 与 A*C 时 间 序 列 进行 
分 析 , 直接 分 解 出 两 序列 的 ~2300 a 周期 (Hallstatt cy- 
cle)IMF 分 量 , 并 进行 互相 关 分 析 , 以 太阳 质心 到 太 
阳 系 质心 的 距离 变化 作为 太阳 轨道 运动 的 指标 , 探 
索 太 阳 轨 道 运动 周期 与 太阳 活动 周期 之 间 的 关系 。 


1 数据 与 方法 


1.1 太阳 质心 到 太阳 系 质心 距离 变化 计算 公式 与 
时 间 序 列 数据 

111 行星 系 质 必 绕 日 会 合 指数 定义 ”经 过 太阳 系 
质心 并 且 垂 直 于 太阳 系 总 角 动 量 矢量 的 平面 为 不 
变 平面 , 它 与 太阳 系 八大 行星 轨道 平面 的 夹 角 平均 
为 1.84。 ,外 行星 的 倾角 都 小 于 1? ,内 行星 除 水 星 倾 
角 较 大 为 6.35$" 外 ,其它 行 星 都 在 2 左右 字 ,因此 可 
以 近似 将 各 行星 轨道 平面 视 为 在 不 变 平面 上 ,由 此 
建立 坐标 系 ,如 图 1 所 示 。 将 太阳 系 中 八大 行 
星 视 为 质点 , P, (i =1,2,3,…,8) 代 表 各 行星 以 质 
量 为 权重 的 位 矢 未 端 , m 、g, 、r, 分 别 为 各 行星 的 
质量 .日 心经 度 以 及 其 相对 于 太阳 质心 的 位 矢 半 
径 。 则 行星 系 质心 (P ) 相 对 于 太阳 质心 (S ) 的 位 
Rr H”: 


rp=SP= = (1) 


注 :S, 太 阳 质 心 ;P, 行 星系 质心 ;P.(i=1,2,3,…,8) 代 表 各 行星 以 质 
量 为 权重 的 位 矢 末端 。 
图 1 行星 会 合 示意 图 


Fig.1 Sketch map of planetary rendezvous 


公式 (1) 为 太阳 质心 $ 到 行星 系 质心 P 的 位 矢 
方程 ,位 矢 的 始 端 为 太阳 质心 ,未 端 为 行星 系 质心 ， 
由 此 可 知 ,行星 系 质心 是 各 行星 以 质量 为 权重 距 太 
阳 质 心 的 权重 半径 的 矢量 和 。 

根据 此 原理 , 刘 复 刚 和 王建 ”把 行星 绕 日 旋转 
轨道 的 权重 半径 矢量 和 的 模 , 定 义 为 行星 会 合 指 
数 。 太 阳 系 行星 会 合 指数 ( KR ) 的 运动 学 方程 为 ”: 


K= 应 sin(g, + o) + Sr cos(g, + oj (2) 


式 中 : K 为 行星 会 合 指数 ; ” 为 行星 轨道 的 权重 半 
径 ; 9, 为 行星 的 初始 日 心经 度 ; o, 为 平均 角速度 ; t 
为 时 间 。 行星 会 合 指 数 可 以 指 代行 星相 会 合 和 相 
背离 的 程度 ,又 因为 行星 以 质量 为 权重 的 半径 x, 的 
矢量 和 代表 行星 系统 质心 的 位 置 , 所 以 行星 会 合 指 
数 玉 还 表示 行星 系统 质心 与 太阳 质心 之 间 的 距离 。 
1.1.2 太阳 质心 到 太阳 系 质心 距离 变化 计算 ”如 果 
把 太阳 系 八大 行星 看 作 一 个 系统 ,根据 二 体 围绕 共 
同 质心 相互 绕 转 的 规律 ,太阳 质心 和 行星 系统 质心 
围绕 太阳 系 质心 相互 绕 转 ,3 个 质心 在 一 条 直线 
上 。 假 设 太阳 和 行星 系统 的 质量 分 别 为 WMs FM, ， 
太阳 质心 和 行星 系统 质心 距离 及 太阳 系 质心 距离 
分 别 为 R M R, o W MsRs=MbR。, 又 由 于 
K=R,+R, ,可 推出 : 
MsRs=M,(K-R,) (3) 
R, = KM,/(M,+M,) (4) 
根据 公式 (2)、(4) ,利用 表 1 中 的 相关 数据 ( 八 
大 行星 的 平均 角速度 ,质量 权重 半径 ,2000 年 1 月 1 
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表 1 太阳 系 行星 轨道 相关 天 文 数据 (历年 =J2000=2000 January 1.5) 


Tab.1 Astronomical data of the planetary orbits in the solar system 


名 称 ”与 不 变 平面 倾角 /” 平均 角速度 had.a: 平均 轨道 半径 /AU ”质量 /10*ks ”质量 权重 半径 /AU ”日 心经 度 /" ”转换 成 弧度 /rad 
水 星 6.35 26.088408 0.38709893 0.33 0.000048 252.25084 4.402608 
金星 25 10.213444 0.72333199 4.87 0.001321 181.97973 3.176145 
地 球 LST 6.283185 1.0000001 1 5.97 0.002238 100.46435 1.753434 
火星 1.63 3.340667 1.52366231 0.64 0.000366 355.45332 6.203831 
AR 0.32 0.529700 5.20336301 1898.19 3.702624 34.40438 0.600470 
土星 0.93 0.213303 9.53707032 568.34 2.031931 49.94432 0.871693 
天 王 星 0.99 0.074790 19.19126393 86.81 0.624538 313.23218 5.466933 
海王 星 0.74 0.038129 30.06896348 102.41 1.154374 304.88003 5.321160 

注 : 第 二 列 和 第 三 列 数据 分 别 引 自 文献 [22] 和 [24]; 其 它 基 础 数据 引 自 NASA。 

日 12:00 的 日 心经 度 , 以 及 其 与 太阳 的 质量 ) 进 行 计 ”edu/,0~13.9 ka BP 数据 分 辩 率 为 5 a,13.9~25.0 ka 


算 , 即 可 获得 太阳 质心 到 太阳 系 质心 距离 随时 间 变 
化 的 数据 序列 。 
1.2 太阳 活动 替代 性 指标 ”C 数据 来 源 

自然 界 中 的 ”C 主要 是 宇宙 射线 中 的 中 子 与 地 
球 大 气 "N 相互 作用 的 产物 。 由 于 地 球 存在 磁场 , 宇 
宙 射 线 在 到 达 地 球 附近 时 会 发 生 偶 转 ,因此 地 磁场 
的 强度 会 影响 大 气 中 的 宇宙 射线 通 量 。 图 2 中 
“C 数 据 的 长 期 趋势 , 正 是 由 于 地 球 磁 偶 极 场 的 变化 
所 导致 的 ,"C 中 的 短期 变动 反映 的 则 是 太阳 活动 的 
变化 “。 太 阳 活 动 强烈 ,太阳 黑子 增加 时 ,会 伴随 
着 更 强 的 太阳 风 。 太 阳 风 的 主要 成 分 是 质子 和 电 
子 ,是 一 种 等 离子 体 囊 电 粒 子 流 ,其 磁 屏 珊 特 性 会 
对 宇宙 射线 产生 较 强 的 作用 ,从 而 使 地 球 大 气 中 的 
“C 生 产 率 降低 。 因 此 太阳 活动 增强 (黑子 数 增多 ) 
对 应 的 是 较 低 的 大 气 “C 生 产 率 ,太阳 活动 减弱 ( 黑 
子 数 减 少 ) 对 应 的 是 较 高 的 大 气 “C 生 产 率 , 这 些 都 


BP 为 10 a,25.0~50.0 ka BP 为 20 a。A"C 时 间 序 列 
中 , 早 于 25.0 ka BP 的 数据 曲线 更 为 平滑 , 旦 误差 相 
对 较 大 ,因此 本 文 研究 集中 于 0~25.0 ka BP。 需 要 
引起 重视 的 一 点 是 ,年 表 中 0~13.9 ka BP 数据 主要 
来 自 树 轮 , 而 更 早 的 数据 主要 来 自 海洋 记录 :珊瑚 、 
有 和 孔 虫 等 。 两 者 的 不 同 在 于 ,树木 可 以 直接 吸收 空 
气 中 的 碳 , 而 珊瑚 等 海洋 生物 吸收 的 溶解 在 水 里 的 
碳 可 能 已 经 被 稀释 相当 长 的 时 间 - > 。 
1.3 EEMD 分 析 方 法 

1998 年 ,Huang 等 ”提出 了 适用 于 非 线 性 、 非 平 
稳 序列 信号 分 析 的 经 验 模 态 分 解 方 法 (Empirical 
mode decomposition, EMD) ,认为 任何 数据 信号 都 是 
一 系列 有 限 的 本 征 模 态 函数 (IMFs) 的 混 乔 。EMD 
方法 ,按照 高 频 至 低频 的 顺序 ,剥离 出 这 些 IMF 分 
量 , 分 解 过 程 不 同 于 小 波 分 析 需 要 借助 基 汪 数 ,而 


被 记录 在 生物 有 机 体 中 ( 树 轮 .珊瑚 等 )。"”C 浓度 随 
时 间 呈 周期 性 变化 ,最 著名 的 是 ~210 a(de Vries cy- 
cle) 和 ~2300 a 周期 (Hallstatt cycle)。 

图 2 显示 50.0 ka BP 的 A"C 记 录 (INTCAL13™)， 
数据 来 源 网 站 http://radiocarbon.webhost.uits.arizona. 


700 25.0 ka BP 


0 10000 20000 30000 40000 50000 
时 间 /a BP 


图 2 大 气 放射 性 碳水 平 A"C 
Fig.2 Radiocarbon series ANC 


是 按照 信号 自身 的 特性 自 适 应 地 把 复杂 数据 序列 
分 解 成 一 系列 具备 平稳 信号 特征 的 IMFs'”。 这 样 
分 解 出 的 每 个 分 量 所 蕴含 的 周期 及 其 物理 意义 便 
更 加 直观 。EMD 方 法 在 广泛 应 用 时 ,发 现 可 能 会 存 
在 模 态 混 徐 ,导致 分 解 出 的 IMF 分 量 可 能 没有 实在 
的 物理 意义 。 针 对 这 一 问题 ,2009 年 Wu 等 ”在 
EMD 基础 上 加 以 改进 ,提出 了 聚合 经 验 模 态 分 解 方 
法 (EEMD ) ,通过 向 信号 中 加 入 高 斯 白 噪 声 来 避免 
模 态 混 释 效应 ,使 最 终 分 解 的 IMFs 保 持 了 物理 上 的 
唯一 性 。EEMD 方 法 能 将 时 间 序 列 分 解 成 若干 真实 
的 有 物理 意义 的 IMF 分 量 和 表示 序列 总 体 趋 势 的 趋 
势 项 

以 EEMD 方 法 直接 剥离 出 太阳 质心 到 太阳 系 质 
心 距离 和 A*C 时 间 序 列 的 Hallstatt 周期 分 量 ,而 不 
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需要 预先 设 定 函 数 关系 式 对 数据 序列 进行 拟 合 , 减 
少 注入 人 为 因素 ,可 以 更 大 程度 上 保留 周期 信号 ， 
提高 计算 精度 。 


2 结果 与 分 析 


2.1 太阳 质心 到 太阳 系 质心 距离 和 A"“C 时 间 序 列 
的 EEMD 分 析 

由 公式 (2) (4) 计 算得 到 0~25.0 ka BP 太 阳 质 
心 到 太阳 系 质心 距离 随时 间 变 化 的 数据 序列 ,数据 
间隔 为 1 a, 绘 制 成 图 3a (为 方便 与 A"*C 数据 作对 
比 ,以 1950 AD 作 为 起 始点 ), 纵 坐标 表示 太阳 质心 
到 太阳 系 质 心 的 距离 大 小 ,单位 为 太阳 半径 ,太阳 
质心 距 太阳 系 质 心 变 化 范围 在 0~2.2 之 间 , 即 最 小 
时 太阳 质心 变化 可 以 与 太阳 系 质心 接近 重合 ,最 大 
时 距 太 阳 系 质心 超过 一 个 太阳 直径 大 小 。 使 用 
EEMD 方 法 对 太阳 质心 到 太阳 系 质心 距离 时 间 序 列 
进行 多 时 间 尺 度 分 解 ,NE=100,Nstd=0.2, 即 集合 样 
本 数 取 值 100, 白 噪声 标准 差 取 原 数 据 序列 标准 差 的 
02 倍 。 分 解 出 高 频 至 低频 的 13 个 IMF 分 量 以 及 一 个 
趋势 项 Res。 各 IMF 分 量 对 应 周期 列 于 表 2。IMF10 
准 周期 为 2273 a, 为 Hallstatt 周 期 分 量 ( 图 3a)。 

使 用 EEMD 方 法 对 0~25.0 ka BP 的 A"C 时 间 序 
列 进行 多 时 间 尺 度 分 解 , 数 据 间 隔 为 10 a, 取 NE= 
100, Nstd=0.2。 分 解 出 高 频 至 低频 的 10 个 IMF 分 量 


以 及 一 个 趋势 项 Res。 各 IMF 分 量 对 应 周期 列 于 表 
2, 计 算出 的 周期 与 文献 ”计算 的 A*C 周期 基本 对 
应 。IMF7 准 周 期 为 2456 a, 为 Hallstatt 周期 分 量 
(图 3b)。 

太阳 质心 到 太阳 系 质心 距离 和 A"*C 时间 序 列 都 
包含 ~2300 a 的 Hallstatt 周 期 分 量 , 将 两 者 绘制 在 图 
3c 中 ,对 比 发 现 ~2300 a 的 Hallstatt 周 期 波动 基本 一 
一 对 应 。 为 进一步 探究 两 序列 Hallstatt 周期 分 量 间 
的 相关 关系 ,利用 IBM SPSS Statistics 23 数据 分 析 软 
件 对 其 进行 互相 关 分 析 。 
2.2 Hallstatt 周期 分 量 间 的 互相 关 关 系 

A“C 年 表 中 0~13.9 ka BP 数据 主要 来 自 树 轮 ， 
而 更 早 的 数据 主要 来 自 海洋 记录 ,因此 以 13.9 ka 
BP 为 分 界线 分 别 计算 0~13.9 ka BP 与 13.9~25.0 ka 
BP 太阳 质心 到 太阳 系 质心 距离 和 A"C 时 间 序 列 
Hallstatt 周期 分 量 间 的 互相 关 关 系 , 如 图 4 所 示 。 

图 4 中 横 坐 标的 正 值 表示 太阳 质心 到 太阳 系 质 
心 距离 变化 领先 于 A"C 变化 , 负 值 则 表示 滞后 于 
A"C。 图 4a 相关 系数 在 相对 位 移 为 320 a 时 达到 极 
大 值 0.52( 在 0.01 级 别 相 关 性 显著 ) ,表示 0~13.9 ka 
BP 太阳 质心 到 太阳 系 质心 距 离 和 A"*C 时 间 序 列 的 
Hallstatt 周期 分 量 间 ,具有 正 相 关 关 系 ,太阳 质心 到 
太阳 系 质心 距离 变化 领先 于 AsC 变化 320 a。 这 与 
Scafetta 等 "中 计算 的 全 新 世 行 星系 质心 相对 于 太阳 
运动 轨道 的 瞬时 偏心 矢量 变化 与 A"C 变化 x/2 (4 


(a) 太阳 质心 到 太阳 系 质心 距离 时 间 序 列 ( 蓝 线 ) 与 其 分 解 出 的 IMF10 分 量 ( 黑 线 ) 
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图 3 太阳 质心 到 太阳 系 质 心 距离 和 AYC 时 间 序 列 的 Hallstatt 周 期 分 量 对 比 图 


Fig. 3 Components with Hallstatt cycle decomposed from the time series of AC and distance between 


the solar centroid and the solar system’s centroid 
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表 2 太阳 质心 到 太阳 系 质心 距离 和 A"C 
时 间 序 列 各 分 量 准 周期 
Tab. 2 Period of the components decomposed from the 
time series of A“C and distance between the solar centroid 


and the solar system’ s centroid 


模 态 分 量 质心 距离 准 周期 /a AMC 准 周期 /a 
IMF1 3.9~12.8 22.6~41.0 
IMF2 12.8 74.0 
IMF3 19.8 130.0~208.0 
IMF4 35.3~45.4 357.0 
IMF5 60.9~92.6 446.0~961.0 
IMF6 171.0~185.0 926.0~2084.0 
IMF7 568.0 2456.0 
IMF8 735.0~1191.0 6253.0 
IMF9 1136.0 25000.0 
IMF10 2273.0 25000.0 
IMF11 6250.0 = 
IMF12 12500.0 一 
IMF13 25000.0 一 

注 : 熏 "表示 无 结果 。 

580 a) 的 相位 差 相 似 。 


图 4b 相关 系数 在 相对 位 移 为 1870 a 时 达到 极 
大 值 0.44( 在 0.01 级 别 相关 性 显著 ) ,这 可 能 是 由 于 
A“C 年 表 中 13.9~25.0 ka BP 数据 主要 来 自 海洋 记 
录 , 不 同 于 树木 可 以 直接 吸收 空气 中 的 碳 , 珊 瑚 等 
海洋 生物 吸收 的 溶解 在 水 里 的 碳 已 经 被 稀释 相当 
长 的 时 间 ' 归 。 因 此 13.9~25.0 ka BP 太阳 质心 到 太 
阳 系 质心 距离 和 A"C 时 间 序 列 Hallstatt 周期 分 量 
间 , 仍 具有 正 相关 关 系 , 但 是 相位 差 有 所 不 同 。 

图 4c 中 0~25.0 ka BP 太阳 质心 到 太阳 系 质心 


(a) 0~13.9 ka BP 
0.8 


(320, 0.52") 


-800 -400 0 400 800 -2300 -1150 


(b) 13.9~25.0 ka BP 
0.8 


距离 和 A"C 时 间 序 列 Hallstatt 周期 分 量 间 的 互相 关 
系数 在 相对 位 移 为 120 a 时 达到 极 大 值 0.35( 在 0.01 
级 别 相关 性 显著 ) ,呈正 相关 。 

以 上 分 析 显 示 , 太 阳 质 心 到 太阳 系 质心 距离 变 
化 与 A"C 变化 均 存在 的 ~2300 a Hallstatt 周期 分 量 呈 
正 相 关 关 系 , 太 阳 质 心 到 太阳 系 质心 距离 变化 可 能 
SEF ANC 变化 320 a。 即 太阳 质心 到 太阳 系 质心 
的 距离 增 大 时 ,大 气 “C 含量 增加 ;太阳 质心 到 太阳 
系 质心 的 距离 减 小 时 ,大 气 剩余 A"C 含量 减 小 。 由 
于 大 气 “C 生 产 率 与 太阳 活动 强度 呈 负 相关 关系 ， 
此 ~2300 a 周 期 中 太阳 质心 远离 太阳 系 质心 时 对 应 
于 太阳 活动 较 弱 时 期 ,太阳 质心 靠近 太阳 系 质心 时 
对 应 于 太阳 活动 较 强 时 期 。 太 阳 轨 道 运动 通过 影 
响 太 阳 活 动 可 能 也 影响 了 地 球 气候 变化 ,如 全 新 世 
0.4~0.6 ka,2.8 ka,5.9 ka,8.2 ka,10.3 ka 的 5 次 冷 事 
件 玉 歼 基 本 发 生 在 太阳 质心 远离 太阳 系 质心 ,太阳 
活动 较 弱 时 期 。 

太阳 绕 太 阳 系 质心 运动 影响 太阳 活动 的 具体 
机 制 尚 不 清楚 ,Scafetta 等 ”认为 太阳 轨道 运动 可 能 
对 太阳 系 动力 学 相关 机 制 具有 指 代 作 用 ,原则 上 ， 
行星 运动 可 以 直接 产生 引力 或 电磁 力作 用 于 太阳 ， 
并 与 日 球 层 相互 作用 ,从 而 调节 地 月 系统 周围 的 宇 
宙 射 线 通 量 以 及 行星 际 宇 家 人 尘埃 浓度 ”3。 刘 复 
刚 "“ 认 为 太阳 自转 角 动 量 与 太阳 绕 太 阳 系 质心 的 
绕 转角 动量 之 和 守恒。 太阳 质心 远离 太阳 系 质心 
时 ,太阳 绕 转 半径 增加 ,引起 太阳 绕 转 角 动 量 增 大 ， 
由 于 太阳 自转 角 动 量 与 绕 转角 动量 之 和 守恒 ,太阳 
自转 角 动 量 减 小 ,导致 太 阳 自 转速 度 减 小 ,太阳 活 


(c) 0~25.0 ka BP 


0.8 


(1870, 0.44”) (120, 0.35") 


0.4 


1150 2300 -800 -400 0 400 800 


时 间 /a 时 间 /a 


注 : 表示 在 0.01 级 别 ( 双 尾 ) 相 关 性 显著 。 


图 4 太阳 质心 到 太阳 系 质心 距离 和 AMC 时 间 序 列 的 Hallstatt 周 期 分 量 互相 关系 数 


Fig. 4 Cross-correlation coefficient of the components with Hallstatt cycle decomposed from the time series 


of A“C and distance between the solar centroid and the solar system’ s centroid 
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动 减 弱 ; 反 之 ,太阳 质心 靠近 太阳 系 质心 时 ,太阳 绕 
转 半径 缩小 ,引起 太阳 绕 转 角 动 量 减 小 ,自转 角 动 
量 增 加 ,从 而 导致 太阳 上 自转 速度 增 大 ,太阳 活动 增 
强 。 上 述 理论 对 于 理解 太阳 质心 到 太阳 系 质心 距 
离 与 太阳 活动 强度 的 负 相 关 关 系 具 有 很 好 的 局 示 
作用 ,但 具体 机 制 仍 要 作 进一步 的 研究 。 


3 结 论 

(1) 文章 基于 行星 会 合 指数 ,推算 太阳 质心 到 
太阳 系 质心 距离 变化 的 公式 ,并 计算 出 时 间 序 列 数 
据 , 发 现 其 变化 范围 在 0~2.2 之 间 ( 单 位 为 太阳 半 
径 ) , 即 最 小 时 太阳 质心 变化 可 以 与 太阳 系 质 心 接 
MEA IANA A DD at —TP ABER 
小 。 利 用 EEMD 方 法 对 0~25.0 ka BP 太阳 质心 到 太 
阳 系 质心 的 距离 时 间 序 列 进行 分 析 , 分 解 出 ~2300 a 
周期 (Hallstatt cycle ) 的 IMF 分 量 。 

(2) 不 同 于 传统 方式 对 于 A*C 时间 序列 的 去 趋 
势 .曲线 拟 合 数据 处 理 ,采用 EEMD 方法 按照 信号 自 
身 的 特性 自 适 应 地 分 解 出 具备 平稳 信号 特征 的 
Hallstatt 周期 分 量 ,使 其 所 蕴含 的 周期 及 物理 意义 
更 加 直观 准确 。 

(3) 由 于 A*C 年 表 中 数据 来 源 不 同 , DA 13.9 ka 
BP 为 分 界线 分 别 计算 0~13.9 ka BP 与 13.9~25.0 ka 
BP 太阳 质心 到 太阳 系 质心 距离 和 A"C 时 间 序 列 
Hallstatt 周期 分 量 间 的 互相 关 关系 ,互相 关系 数 可 
达 0.52、0.44, 呈 正 相 关 。 太 阳 轨 道 运 动 与 太阳 活动 
的 关系 研究 不 再 局 限于 全 新 世 , 拓 展 到 25.0 ka BP. 

(4) 以 太阳 质心 到 太阳 系 质心 距离 变化 指 代 太 
阳 轨 道 运 动 ,以 A"*C 变化 指 代 太 阳 活 动 , 两 者 均 存 在 
~2300 a 的 Hallstatt 周 期 分 量 , 并 有 旦 具有 相关 性 ( 太 
阳 质 心 到 太阳 系 质心 距 离 变 化 与 A"C 数 据 呈 正 相 
关 , 即 与 太阳 活动 变化 星人 负 相 关 ), 表 明太 阳 质 心 远 
离 太 阳 系 质心 时 ,大 气 剩 余 A"C 含 量 增加 , 即 太 阳 活 
动 减弱 ;太阳 质心 靠近 太阳 系 质心 时 ,大 气 剩余 A"C 
含量 减 小 , 即 太 阳 活 动 增强 。 
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Abstract: Solar activity, evidenced by “C proxies, shows a period of ~2300 years (Hallstatt cycle). Its physical 
origin remains uncertain. Recent studies suggested that ~2300-year Hallstatt oscillations of solar activity may be 
caused by the solar inertial motion (the Sun motion around the center of mass of the solar system). In this study, a 
formula for calculating the distance between the solar centroid and the solar system’ s centroid is deduced based 
on the kinematic equation of the planet juncture index. In the appropriate proxy of solar inertial motion, the time 
series of distance change is reconstructed using this formula. Ensemble empirical mode decomposition (EEMD) 
is a method with which any complicated signal can be decomposed into several intrinsic mode functions (IMFs). 
This decomposition method is adaptive and highly efficient. It applies to nonlinear and non-stationary processes 
since the decomposition is based on the local characteristic time scale of the data. To compare the radiocarbon 
series of INTCAL13 data with the variations of the distance between the solar centroid and the solar system’ s 
centroid from 0 to 25.0 ka BP, we employ the empirical mode decomposition method. The analysis results are as 
follows: (1) Various components have been identified by decomposing the distance signal, including a significant 
component with a period of ~2300 years. (2) Following the analysis of the time series of distance change, there is 
also a component with ~2300-year Hallstatt cycle in the data on the production rate of radiocarbon. The traditional 
method of A“C record preprocessing, applying a pre-selected detrending function, amounts to injecting external 
information into data. Different from the traditional method the true and more physical meaningful component 
with the Hallstatt cycle is decomposed by the adaptive EEMD method. (3) Next, the ~2300-year variations of the 
two series considered are compared. The last 13900 years data in INTCAL13 record derive mostly from tree-ring 
chronologies, while the older data is made using mostly marine records. Thus, cross-correlation analysis is applied 
to the above two time-frames of the components with the Hallstatt cycle decomposed from the time series of the 
distance and A“C. The cross-correlation coefficient is 0.52 from 0 to 13.9 ka BP, and 0.44 from 13.9 to 25.0 ka 
BP. The values confirm the assumption that the variations of the distance between the solar centroid and the solar 
system’ s centroid and the production rate of radiocarbon are related. As discussed above, the study of the 
relationship between solar inertial motion and solar activity covers the Holocene and extends to 25.0 ka BP. (4) It 
appears that solar inertial motion is the cause of the Hallstatt cycle found in solar activity proxy. The “C 
concentration of the atmosphere increases when the solar centroid moves away from the solar system’ s centroid, 
which corresponds to the period of weak solar activity. When the solar centroid is near the solar system’ s 
centroid, the “C concentration decreases which corresponds to the period of strong solar activity. 


Key words: solar inertial motion; solar activity; Hallstatt cycle; ensemble empirical mode decomposition 
(EEMD) 


